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Abstract 

We investigate the exact relation existing between the stability equation for the solutions of a 
mechanical system and the geodesic deviation equation of the associated geodesic problem in the 
Jacobi metric constructed via the Maupertuis-Jacobi Principle. We conclude that the dynamical 
and geometrical approaches to the stability /instability problem are not equivalent. 

1 Introduction 

In recent years, several authors [1], [2], [3], [4], have formulated geometrical criteria of (local) stabil- 
ity/instability in mechanical systems using different "geometrization" techniques (Maupertuis-Jacobi 
Principle, Eisenhart metric, etc). The main idea is to interpret the local instability problem, under- 
stood in terms of sensitive dependence on initial conditions, as the study of an appropriate geodesic 
deviation equation. As a principal application, chaotic behaviors in Hamiltonian mechanical systems 
that appears in cosmological models have been described using these results. Most of these works 
are constructed using the Maupertuis-Jacobi principle for natural mechanical systems, both in the 
very well known Riemannian case, but also in the recent generalization to the non-Riemanian one 
[2]. 

The Maupertuis-Jacobi principle establishes, in his classical formulation, the equivalence between 
the resolution of the Euler-Lagrange equations of a natural Hamiltonian dynamical system (hence 
the Newton equations), for a given value of the mechanical energy, and the calculation of the geodesic 
curves in an associated Riemannian manifold. Throughout the time this equivalence has been used 
for different purposes, as the mentioned description of chaotic situations, but also in the analysis 
of ergodic systems [3], [5], non-integrability problems [7], determination of stability properties of 
solitons [8], [9], etcetera. 

The linealization of the geodesic equations in a given manifold gives in a natural way the so-called 
Jacobi equation, or geodesic deviation equation, that allows to compute the stability /instability of a 
given geodesic curve in terms of the sign of the curvature tensor over the geodesic (in fact, for two- 
dimensional manifolds, the problem reduces simply to the computation of the sign of the gaussian 
curvature along the geodesic, see for instance [6]). 
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The geometrization of the mechanical problem provides, as mentioned, a possible criterion of 
stability of the solutions of in terms of the geodesic deviation equation of the Jacobi metric associated 
to the system, via the Maupertuis-Jacobi principle, that we will call Jacobi-metric stability criterion. 

From the point of view of the Variational Calculus applied to geodesies, a similar result is obtained 
for the problem of calculation of fixed-endpoints geodesies, where the sign of the second variation 
functional is determined by the geodesic deviation operator. 

In this work we analyze the exact relation existing between this Jacobi-metric criterion and the 
direct analysis of the stability of the solutions without using the geometrization principle. The 
linealization of the Euler-Lagrange equation (in this case, Newton equations) lead to a Jacobi-like 
equation that generalizes the geodesic deviation one to the case of natural mechanical systems. In 
fact, this equation is also called Jacobi equation in the context of second-order ordinary differential 
equations theory or KCC theory (Kosambi-Cartan-Chern), [10], [11]. 

As we will see, the two approaches (geometrical and dynamical) are not equivalent in general, and 
the Jacobi-metric criterion do not provide exactly the same result as the standard (or dynamical) 
one. 

The structure of the paper is as follows: in section 2 we present the concepts involved in the 
work; Section 3 is dedicated to Jacobi-metric stability criterion and its relation with the dynamical 
one. In Section 4, the analysis is extended to the variational point of view for fixed end-points 
problems. Finally, an Appendix is included with several technical formulas (more or less well known) 
about the behavior of covariant derivatives and curvature tensor under conformal transformations 
and reparametrizations of curves. 

2 Preliminaries and Notation 

We treat in this work with natural Hamiltonian dynamical systems, i.e., the triple (M,g,C), where 
(M, g) is a Riemannian manifold, and C is a natural Lagrangian function: C : TM — > R, £ = T — U, 



in a system of local coordinates (q 1 ,...,q n ) in M, U is a given smooth function U : M — > R, 
= (o ,1 (i), • • • , Q n {t)) is a smooth curve on M, and gij are the components of the metric g in this 
coordinate system (Einstein convention about sum in repeated indices will be used along the paper). 

The solutions (trajectories) of the system are the extremals of the action functional S[j], defined 
in the space of smooth curves on M: 7 : [to,h] — > M, (we assume that 7 is at least of class C 2 in 
the interval (to,ii)). 



T = 





to 



(1) 
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where 7 € T(TM) stands for the tangent vectorfield i.e. j(t) — ^-(t) G T^n\M. 

Euler-Lagrange equations associated to this functional are Newton equations for the system: 

SS = =► V^7 = -graded (2) 

where V-y stands for the covariant derivative along j(t) = (q l (t)): 

being T l - k the Christofell symbols of the Levi-Civitta connection associated to the metric g. 

r i = 1 ki f ^9jk d(hk _ dgij_ 
ij 2 9 V dq i dqj dq k 

gradC/ is the vectorfield with components: (gradC/)* = fl^^j- Equation (2) is thus written in local 
coordinates as the following system of ordinary differential equations: 

f=t+Wt = -j>% (3) 

Natural Hamiltonian dynamical systems over Riemannian manifolds satisfy Legendre's condition 
in an obvious way, and thus the Legendre transformation is regular, i.e. there exists a diffeomorphism 
between the tangent and cotangent bundles of M in such a way that the Euler-Lagrange equations 
are equivalent to the Hamilton (or canonical) equations. 

oh . oh 

P * = ~W q (4) 

where 

Pj = W = 9ij ^'' H = 2 9l3piPj + U 
and g %3 denotes the components of the inverse of g. 

This kind of systems are autonomous, thus the mechanical energy is a first integral of the system: 

E= l -g ij q i tf+U{q\...,q n ) 

Stability of the solutions of (3), understood in terms of sensitive dependence on initial conditions, 
is interpreted as follows: The trajectory j(t), solution of (3), is said to be stable if all trajectories 
with sufficiently close initial conditions at to remains close to the trajectory j(t) for later times t > to- 

Let j(t; a) = a), . . . , q n (t; a)) be a family of solutions of equations (3), with j(t) = j(t; 0), 
and given initial conditions q % (to; a), q l (to;a). Let us assume that the initial conditions are analytic 
in the parameter a. Then: j(t) = (q l (t)) is a stable trajectory if for any e > 0, there exists a 5(e) > 
such that \q l (t;a) — q l (t)\ < e for t > to and for all trajectories q(t;a) = (q l (t;a)) satisfying both 
|^(t ;a) - (/'(to)! < S and \q* (t ; a) - c? (t )| < 5. 
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Assuming that g is smooth and considering that j(t; a) are analytic in a (they are solutions of 
an analytic system of differential equations), we can write, for a sufficiently small: 



q%a) =q i (a) + av i (t) + o(a 2 ) , 



da 

In a similar way, we can write: 



(5) 



a=0 



<9P 

Y) k {q{t-a)) = T) k (q{t)) + a-^(q{t))v\t) + o(a 2 ) (6) 



g ij (q(t;a)) = (q(t)) + a^(q(t)) v l (t) + o(a 2 ) (7) 

djU(q(t;a)) = djU{q{t)) + a didjU(q(t)) v l (t) + o(a 2 ) (8) 

where djU = 0. 

Thus equations (3) become: 

v* + 2T) k viq k = -g^v 1 (d^U + Y\ p djU) + g^djUv 1 (9) 

where all functions are taken at Taking into account the expression of the second order covariant 
derivatives: 

^~=v l + d l Y) k q l q^v k + 2T) k v^ q k + Y' k .c'if + T\ p T l ]k q^cfv k 

and the components of the Riemann curvature tensor: R(X,Y)Z = — Vx(VyZ) + Vy(Vx-Z') + 
V [XX] Z, VX,Y,Z £T(TM): 



we finally arrive to the expression: 

D 2 v i 



+ Rj jk q l q j v k = -g* {d^U - T^U) v l 



dt 2 

that can be written as a vector equation: 

VjVjV + KA j {V) + Vygradf/ = (10) 
where V = V(t) = (v l (t)), and we have used the sectional curvature tensor: 

K x (Y) = R(X, Y)X, VX, Y G F(TM) 
and the Hessian of the potential energy U: 7~t(U) = VdU 

VdU = (djdiU - d k UY k ^j dq j ® dq l 
in such a way that VX, Y € T{TM) 

VdU(X,Y) = (Vxgrad(C/),Y) = (Vygrad(f/), X) 



Solutions of equation (10) determine the behavior of the family of solutions j(t, a) with respect 
to the selected solution ^y(t). Thus typical solutions of linear equations (trigonometric functions, 
exponentials, etc.) will prescribe the stability/instability situations. In several contexts equation 
(10) is usually called Jacobi equation, by analogy with the geodesic case. In fact, in the so-called 
KCC theory on second order differential equations, equation (11) is nothing but the Jacobi equation 
for the special case of Newton differential equations. In order to avoid confusions we will denote 
Hessian operator for the mechanical system to: 

AV = VjVjV + K^V) + Vygi&dU 

and thus we reserve the term Jacobi operator (and equation) to the geodesic case, i.e. to the geodesic 
deviation equation. 

In the special case of fixed starting point for the family of solutions j(t;a), i.e. j(to;a) = 
7(to), an equivalent approach to equation (11) can be considered. The first variational derivative 
of functional (1) lead to Euler-Lagrange equations (3), and thus the second variation functional 
(or Hessian functional) will determine (together obviously with the Legendre straightness condition, 
automatically satisfied for this kind of systems, see [12]) the local minimum/maximum character of 
a solution of (3). The second- variation functional of the action S, for the case of proper variations 
(V G Y(TM) such that V(t ) = V(ti) = 0) is: 

S 2 S[j{t)} = - / 1 dt (V^V^ V + K^V) + VygradC/, V) = - I ' dt (AV, V) (11) 

J to J to 

and thus the positive or negative definiteness of the A operator determines the character of the 
solution j(t). 

3 The Jacobi-Metric Stability Criterion 

The Maupertuis-Jacobi Principle establishes the equivalence between the resolution of the Newton 
equations (3) of the natural system and the calculation of the geodesic curves in an associated 
Riemannian manifold. The crucial point of the Principle is the existence of the mechanical energy as 
first integral for equations (3). Solutions of (3) corresponding to a fixed value E = T + U will be in 
one to one correspondence with the solutions of the equations of geodesies in the manifold M with 
the so-called Jacobi metric: h = 2(E — U)g, associated to the E value. 

Geodesies in the Riemannian manifold M = (M, h) 1 can be viewed as extremals of the free-action 

1 We will call s = Sh, i.e.: ds 2 g — gijdq l dq\ ds 2 — h i jdq 1 dq j , and hij = 2(E — U)gij. We will also write V J 
for the covariant derivative with respect to h, and, for any vectorfields X, Y G T(TM): h(X,Y) = {X,Y) J , and 

\\x\\ j = ^/(x,xy. 
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functional Sq or of the Length functional L: 

Sob] = ^(||7(i)H J ) 2 ^; m = r ll7(*)H J ^ (12) 

Jto z Jt 

for any differentiable curve 7 : [to, h] — > M connecting the points 7(^0) = P an d t(*i) = Q, P,Q £ M. 
The extremal conditions, SSq = and 5L = 0, lead us to the Euler-Lagrange equations (equations 
of the geodesies in M): 

5S = => V^ 7 = 0; 5L = V( 7 = A(t)<y, A(t) = ^-J (13) 

= leads to the equations of the geodesies parametrized with respect to an arbitrary parameter 
t (often called pre-geodesics) as a natural consequence of the invariance under reparametrizations of 
the Length functional, whereas SSq = produces the equations of affmely parametrized geodesies. 
If we restrict to the arc-length parametrization and we will denote, as usual, 7' = equations 
(13) are written as: VyY = 0, or explicitly, in terms of Christoffel symbols Y l - k of the Levi-Civitta 
connection of h, as: 



D(q 



=( q r + T) k ^y^y = o (14) 

The Maupertuis-Jacobi Principle can be formulated in the following form: 
Theorem of Jacobi. The extremal trajectories of the variational problem associated to the func- 
tional (1) with mechanical energy E, are pre-geodesics of the manifold (M, h), where h is the Jacobi 
metric: h = 2(E - U) g. 

From an analytic point of view, the theorem simply establishes that the Newton equations (3) 
for the action S, are written as the geodesic equations in (M,h): Vy7 ; = 0, when the conformal 
transformation: h = 2{E — U) g, and a reparametrization (from the dynamical time t to the arc-length 
parameter s in (M,h)) are performed. 

Moreover, the dependence between the two parameters is determined over the solutions by the 
equation: 

- = 2^E-U{ 1 (s))T = 2(E - U( 7 (s))) (15) 

The proof of this theorem can be viewed in several references (see for instance [6], see also [12] for 
a general version of the Principle). However, a very simple proof of the theorem can be carried out 
by the explicit calculation of equations (14) in terms of the original metric g, making use of Lemmas 
1 and 2 of the Appendix, that detail the behavior of the covariant derivatives under conformal 
transformations and re-parametrizations. Vy7' = turns out to be 

VyY + (grad(ln(2(£ - [/))), Y)Y " Y>grad(ln(2(£ - 17))) = (16) 
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in terms of the V derivative. By applying now Lemma 2 to (16) we obtain, after the corresponding 
reparametrization and simplifications, the equation 

Vy7 + grad(£/) = 

i.e. the Newton equations of the mechanical system. 

This result allows to define the Jacobi-metric criterion for stability of the mechanical solutions in 
terms of the corresponding geodesies of the Jacobi metric. 

In an analogous way to the previous section, one can linearize the equations (14) of the geodesies 
in (M, h) by considering a family of geodesies 7(5; a): 

7(t;a) = 7(t) + aV + o{a 2 ) 

with V{s) = -^fe-^ • Following the same steps, one finally arrives to the expression 

V$V$V + K$(V) = (17) 

where V = V(s) = (v l (s)), and K J is the sectional curvature tensor of the h metric. 

Equation (17) is the Geodesic Deviation Equation, or Jacobi Equation, for a given geodesic 7(5) 
of (M,h). We will denote Jacobi Operator, or Geodesic Deviation Operator to: 

A J V = Vy V J y V + K$(V) (18) 

Thus stability of a solution of Newton equations j(t) will be determined, in this criterion, if the 
corresponding geodesic 7(5) is stable, that finally leads to equation (17). 

In order to determine the exact relation existing between the Jacobi-metric criterion and the 
dynamical o standard one, we will analyze now equation (17), by using the results about conformal 
transformations and re-parametrizations included in the Appendix. 

Applying Lemma 1 and Lemma 3 (see Appendix) to the Jacobi operator (18) and simplifying the 
expressions, equation (17) is written as: 

A J V = V y V Y V + K Y (V) + ^(F,V)V 1 n' + (F, 1 ')V Y V-^( 1 ', 1 ')V v F + 

+ (<F, Vy V) + ± (F, V) (F, i) + (VyF, 7')) i + 

+ Q (f, vy y> + \ (f, y> 2 - \ W, 7') (F, F>) V + 

+ (~ (VyV, v) - (y, v y y> - \ (f, y> <y, v)\ f (19) 

depending only on the metric g, and where F denotes: F = gradln(2(£' — U)). Re-parametrization 
of 7(5) in terms of the i-parameter: 
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and application of Lemma 2 to (19) lead to: 

A ' ,V = (2(E - U)Y ( V ^ y + K ^ V ) + \ ^ y ) V 77 ~ \ (7. 7) VyF+ 



+ ((F,V^) + (VyF,7))7 + 



+ (-^(^V^ 7 )-(7,V^})f) (20) 

Expression (20) is written in terms of quantities depending only on the metric g and the t- 
parameter. In order to relate this expression with the Hessian operator A we need to remember 
that j(t) is a solution of the Newton equations (3) of energy E, and thus: V^j = — gradf/, (7, 7) = 
2(E — J7( 7 (i))). Using these facts and simplifying we arrive to: 



A J V 



d ((V,gi&dU)\ (graded, V) + (7, V^V) 



AV -dt(^Trn + E-u 1 ^ adu 



(21) 



(2(E-U)y 

where we have used the identity: (7, X/ygr&dU) = (V, V^gradC/). 

Obviously, the two operators do not coincide, and correspondingly solutions of the Jacobi equation 
A J V = and the equation AV = do not so. The two criteria of stability are not equivalent. In 
order to investigate equation (21) to determine the reasons of this non-equivalence between the two 
criteria, we have to remark that whereas all the geodesies 7 (s; a) considered in the calculation of 
correspond to mechanical solutions of energy E (they are solutions of the equation of geodesies in 
(M, h), with h = 2(E — U)g), the solutions j(t; a) are in principle of energy: 

E a = ^{t-^gij^t-a^it-a) + U(q(t;a)) (22) 

But a correct comparison between two stability criteria is only well established if the criteria act 
over the same objects. Thus the comparison is only licit if one restricts the family j(t; a) to verify: 
E a = E. Expanding (22) in a we find: 

E a = E + a « 7 , V^V + (grade/, V})) + o(a 2 ) (23) 

And thus the requirement E a = E reduces to the verification of: (7, V-y) = — (gradf/, V). 
Thus the relation between the Jacobi operator and the hessian operator restricted to equal-energy 
variations is: 



A J V = 



dt\ E-U 1 



(2(E-U)Y 

and the two operators are not equivalent, even considering the equal-energy restriction. 



(24) 
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4 The Variational point of view 



As it has been explained in the Introduction of this work, we will apply now the above obtained 
results to the special case of fixed end-points, i.e. we will restrict our analysis to the situation where 
the conditions: 7(to) = P an d 7(*i) = Q, with P and Q fixed, are imposed . From the mechanical 
point of view, this is exactly the case of the calculation of solitonic solutions in Field Theories (see 
for instance [8]) where asymptotic conditions determine the starting and ending points. Using the 
Maupertuis-Jacobi Principle, this situation is translated to the problem of calculating the geodesies 
connecting two fixed points in the manifold M. We thus use the framework of the Variational 
Calculus for fixed end-points problems. 

The minimizing character (local minimum) of a geodesic 7(5) connecting two fixed points is 
determined by the second variation functional: 

S 2 S = - f 1 (A J V, V) ds , 5 2 L = - r /A J V^, ds (25) 

J so J so 

where A J is the geodesic deviation operator of h: 

A J V = V y V y V + R J (j', V)j' = Vy Vy V + Ky (V) 

where V € T(TM) denotes any proper variation and V ± is the orthogonal component of V to the 
geodesic. 

We will show now two theorems, in the first one it is established the difference between the second 
variation functional of the dynamical problem and the corresponding one to the free-action functional 
associated to the Jacobi metric. In the second one, a similar analysis is carried out for the Length 
functional. 

Theorem 1. Let j(t) be an extremal of the functional 5 [7] = (^ (7,7) — U(j)) dt, and let 
Soil] = Jsq 1 \ (7' ' il') J ds be the free-action functional of the Jacobi metric associated to 5 [7] and 
corresponding to a fixed value, E, of the mechanical energy, then the corresponding Hessian func- 
tionals verify: 

6 2 S{[ 1 ] = 5 2 S[j] + / 1 dt 2 (7, V^V) (F, V) (26) 
J to 

where F = gradln(2(£ - U)). 

Theorem 2. Let ^y{t) be an extremal of the S[j] = (^(7>7) — ^(7)) dt functional and let 
L J [l] = \\l'\\ ds be the length functional of the Jacobi metric associated to £[7] and corresponding 
to a fixed value, E, of the mechanical energy, then the corresponding hessian functionals verify: 

8 2 L J [ 7 ] = 5 2 S[ 7 ] - J 2{E a _ - [(Vyr, V) - s (7, V^V)} 2 (27) 
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From (27) it is obvious that minimizing geodesies are equivalent to minimizing (stable) solutions of 
the dynamical system, i.e. a positive defmiteness of 5 2 L J implies the same behaviour for S 2 S, but it 
is not necessarily true the reciprocal statement. 

If we restrict the variations to the orthogonal ones, V = V ± , (27) can be re- written as: 

6>s\ v=v± = 6*l j + r ds ((^,yM 2 

J so 

The proofs of these two theorems are based on the behaviour of the covariant derivatives and the 
curvature tensor under reparametrizations and conformal transformations of the metric tensor. We 
thus use the technical results included in the Appendix. 

Proof of Theorem 1. We start with equation (25) particularized to the case of the Jacobi metric: 

r>Si 



^oM= ds(-A J V,VY 



with A J V = V^V^V + K$(V). 

Using expression (21), deduced in the previous section after changing the metric and re-parametrizing, 
we can write: 

(y J y VifV + K],(V), V) J = *_ <V^V + K,,(V) + VvgradC/, V) + (28) 

+ Ke L u) I {{F ' ¥){A " v)) + TwhlT^ ViVXvsdu, v) 

And thus, the second variation functional is written as: 



2 (0) = - r ds(V J Y V J y V + K$(V),V) 

J sn 



so 

= - / 1 dt (V^V^V + Ka,(V) + Vygradf/, V) + 

+ f 1 dt2(^,V i V)(F,V)- (F,V) (j,V)\% 
J to 



For proper variations: V{t\) = Vfo) = 

with F = gradLn(2(£ - U)) = graded 

E — U 

Q.E.D. 

Proof of Theorem 2: For the Length functional we have: 



^P(O) = - ds (V^V^ + K^), V±y 



10 



where 
and thus: 

By using Theorem 1, we have that 



^2 



K0)+ [ 2 2(^,V^V)(F,Vdt) + f 2 A(t)dt 



where: 

f 2 A(t)dt = - rds(( 1 f ,V^V) J ) 2 = - f t2 dt(2(E-U)f( 1 ',V J ll V) 2 
From Lemma 1 and Newton equations, we have 

J* 2 A (t)dt = - J* dt 2{E l _ v) ((7, V^V) + {E- U)(F, V)) 2 

Finally 
Q.E.D. 



(o) = ^|P(0) - jT *2(^tf) [<-r> v ^> " ^7, v)Y 



5 Appendix 

Lemma 1. Given a conformal transformation in a riemannian manifold, (M, g) — > (M, g); g = /(x) 5, 
/(x) / 0, Vx € M, let V and V be the associated Levi-Civita connections respectively. Then, for all 
X,Y,Z G r(TM) it is verified that: 

V X Y = V X Y + 1 -{F, Y)X + l -(F, X)Y - l -(X, Y)F (29) 

W X V Y Z = V x V Y Z+^(F,Z)VxY + ^(F,Y)VxZ-^(Y,Z)V x F+^{F,X)V Y Z + 
+ Q(F, V y Z) + 1 -{F,Z)(F,Y) \(Y,Z)(F,F))X + 

+ Q(VxF, Z) + \{F, V X Z) + \{F, Z)(F, X)^jY + 

+ Q(VxF, y) + l(F, V X Y) + \{F, X)(F, F>) Z + (30) 

+ (y(v^,z) - ±<y, VxZ) - i(x,VyZ) - l(F,z)(x,y) - i<F,y><x,z}) f 

where the scalar products are taken with respect to the metric g and F = grad (In /) (grad stands 
for the gradient with respect the metric g). 
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Proof: By direct calculation. Let us consider the expression of the Christoffel symbols of the g metric: 

1 .V 

and substitute g^ = f gtj, g y = jg v ' ■ So 



fjfe = T,~9 ir (d k g jr - d r g jk + djg rk ) 



2 

and the covariant derivative will be 



fjfe = T) k + \ [grad(ln./T {5)g mk + 5\g mo ) - 5jfe grad(ln./T 



v x y = x^y = x> + r} fc y*) ^ = v*y + 4^*^ 



where A*- fe stands for: 



A) k X*Y k = \ [grad(ln/) m (8)g mk + Slg mj ) - «fe fc grad(ln /)*] X*Y k 
= \ «grad(ln/),y>X< + (grad(ln/), X)Y*) - ±(X, r)grad(ln/) J 



Finally, simplifying 

VxY = V X Y + i(grad(ln/),y)X + l(grad(ln/),X)y - l(X, y>grad(ln/) 

and similarly for (30). 
Q.E.D. 

Lemma 2. Given a (differentiable) curve 7 : \t\,t-^ ^ M on M, let 7(5) = j(t(s)) be an admissible 
re-parametrization of 7, ds = /(x(t))di f/(x(t)) / 0,Vt € [ti,t 2 ]j- Tiien € T(TM): 

1 



Vyl = -jy-^^^X (31) 



Vy 7 ' = (V77 ~ (grad(ln /), 7)7) (32) 

VyVyX - -JL (V^X - (grad(ln/), 7 ) Vyf) (33) 



wiiere j(t) = and 7 '(s) = ^f) 
Proof: Again by direct calculation 



r, / lo , 1 /da;' 1 . . , fc \ 9 /Id /i 4 \ .,-. fc l\ 9 

= 72 ( V ^7 - 9 fc In . /•.;•' . K^;) = ^ (Vy7 - (grad(ln/), 7 ) 7 ) 



Vy VyX = Vy (j V; ( Xj = Vy Qj V^X + V^X = 

= £ I ()) + 72 = 7^ (V ^ X - (grad(ln/), 7 ) V^X) 
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Q.E.D. 

Lemma 3. Given a conformal transformation in a Riemannian manifold: (M,g) — > (M,g), g = 
f(x)g, let R and R be the associated curvature tensors respectively. Then, for any X,Y, Z G T(TM), 
it is verified that: 

R{X,Y)Z =R{X,Y)Z -l(X,Z)V Y F+^(Y,Z)WxF + 

+ Q (VyF, Z)-\ (F, Z) (F, Y) + \ (Y, Z) (F, F>) X + 

+ (-1 (V X F, Z) + \ (F, Z) (F, X)-- A (X, Z) (F, F>) Y + 

+ Q (V Y F, X)-± (V X F, Y)^j Z + (34) 

+ (±(F,Y)(X,Z)-±(F,X)<y,Z?) F 

where V is the Levi-Civita connection associated to g, F = grad(ln/) and the scalar products and 
the gradient are taken with respect to the metric g. 

Proof: Apply Lemma 1 to the formula: R(X, Y)Z = -V x (VyZ) + Vy(V x Z) + V [x ,y]Z, and simplify. 
Q.E.D. 
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